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Abstract 

The particle-bed core offers mitigation of some of the problems of solid-core nuclear 
rocket reactors. Dividing the fuel elements into small spherical particles contained in a 
cylindrical bed through which the propellant flows radially, may reduce the thermal 
stress in the fuel elements, allowing higher propellant temperatures to be reached. The 
high temperature regions of the reactor are confined to the interior of cylindrical fuel 
assemblies, so most of the reactor can be relatively cool. This enables the use of 
structural and moderating materials which reduce the minimum critical size and mass of 
the reactor. One of the unresolved questions about this concept is whether the flow 
through the particle-bed will be well behaved, or will be subject to destructive flow 
instabilities. Most of the recent analyses of the stability of the particle-bed reactor have 
been extensions of the approach of Bussard and Delauer, where the bed is essentially 
treated as an array of parallel passages, so that the mass flow is continuous from inlet to 
outlet through any one passage. A more general three dimensional model of the bed is 
adopted here, in which the fluid has mobility in three dimensions. Comparison of results 
of the earlier approach to the present one shows that the former does not accurately 
represent the stability at low Re. The more complete model presented here should be 
capable of meeting this deficiency while accurately representing the effects of the cold 
and hot frits, and of heat conduction and radiation in the particle-bed. It can be extended 
to apply to the cylindrical geometry of particle-bed reactors without difficulty. From 
the exemplary calculations which have been carried out, it can be concluded that a 
particle bed without a cold frit would be subject to instability if operated at the high 
temperatures desired for nuclear rockets, and at power densities below about 4 
megawatts per liter. Since the desired power density is about 40 megawatts per liter, it 
can be concluded that operation at design exit temperature but at reduced power could be 
hazardous for such a reactor. But the calculations also show that an appropriate cold frit 
could very likely cure the instability. Mora definite conclusions must await calculations 
for specific designs. 

* R.C. Maclaurin Professor of Aeronautics and Astronautics 


NTP: Technology 


498 


NP-TIM-92 



Conclusions 


Comparison of three quite different approaches to modeling the stability of the particle-bed 
reactor, all with consistent physical assumptions, shows that a complete linear stability such 
as that presented here is in fact necessary for reliable prediction of the stability of the 
particle-bed reactor. The approach, termed here the Parallel -Stream model, where the 
reactor is assumed to be composed of a series of channels coupled only at their inlets and outlets, 
does not accurately represent the stability at low Re, nor does it represent the effect of heat 
conduction in the bed. 

The model termed here (perhaps somewhat naively) the Complete Model should be capable of 
accurately representing the effects of the cold and hot frits, and of heat conduction and 
radiation in the particle bed. It can be extended to apply to the cylindrical geometry of 
particle-bed reactors without difficulty. 

From the exemplary calculations which have been carried out, it can be concluded that a 
particle bed without a cold frit would be subject to instability if operated at the high 
temperatures desired for nuclear rockets, and at power densities below about 4 megawatts per 
liter. Since the desired power density is about 40 megawatts per liter, it can be concluded that 
operation at design exit temperature but at reduced power could be hazardous for such a reactor. 
But the calculations also show that an appropriate cold frit could very likely cure the 
instability. More definite conclusions must await calculations for specific designs. 
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SCHEMATIC OF PARTICLE-BED 



centerline of annular fuel element 
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GOVERNING EQUATIONS 


Vp-£s - ^1 u*| pu (I) 

P _ |S() 0 


b _ , 7S <] CJ 1 
K" ' r* ”p 

(2) 

3T P , 

Ps <>s ; H P»«p H p h * k cfr v r p 

(3) 

k cff = k coiul ■* k rsui ' k co!i(l 1 a ^p 1 p^ 

H) 

h., to <’ n? l + .875 (l V 

tT p r Z Kc Op r 7- 

to) 

•rr- . 

l K: p ^ 1 P*' ' l n,t p ( *p 1 |',)b 

K») 

♦ V (pu) 11 
ill 

(71 

l» l»«* 

(Hi 


NP-TIM-92 


501 


NTP. Technology 


NON-DIMENSIONAL GOVERNING EQUATIONS 
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ZEROTH-ORDER OR STEADY SOLUTION 
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FIRST ORDER 
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Variables are; p, p, Tg, u and T p 
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PARAMETERS 

Dimensionless Parameters : 

from In, v, b 1( b 2 
from 3n, q. c. H, K c , K r 

Operating Parameters : 

T g0 (exit) - 3000 K 
po (exit) - 100 atm 
Q - 4x1 0 10 watt/m3 

Design Parameters : 

I - o 01 m 
D p - .5 x 10-3 m 

Stability Parameters : 

_ P()(0) up ffi) p p __ PpQI 

c ~ m 
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APPROACHES TO INSTABILITY ANALYSIS 

1) Parallel Flow Instability 

2) Local Instability Analysis 

3) Full Stability Analysis 


PARALLEL-STREAM INSTABILITY 


Instability is possible it po(1) increases with mass flow 
density for fixed Q and po(0). Hence: 
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COMPLETE INSTABILITY MODEL 
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LOCAL INSTABILITY 
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COMPLETE INS I ABILITY MODEL 
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